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ABSTRACT 


The Oort Cloud Comets (OCCs), exemplified by the Great Comet of 1997 
(Hale-Bopp), are occasional visitors from the heatless periphery of the solar sys- 
tem. Previous works hypothesized that a great majority of OCCs must physically 
disrupt after one or two passages through the inner solar system, where strong 
thermal gradients can cause phase transitions or volatile pressure buildup. Here 
we study the fate of small debris particles produced by OCC disruptions to de- 
termine whether the imprints of a hypothetical population of OCC meteoroids 
can be found in the existing meteor radar data. We find that OCC particles 
with diameters D < 10 //m are blown out from the solar system by radiation 
pressure, while those with D > 1 mm have a very low Earth-impact probability. 
The intermediate particle sizes, D ~ 100 /jm, represent a sweet spot. About 
1% of these particles orbitally evolve by Poynting-Robertson drag to reach or- 
bits with semimajor axis a ~ 1 AU. They are expected to produce meteors with 
radiants near the apex of the Earth’s orbital motion. We find that the model 
distributions of their impact speeds and orbits provide a good match to radar 
observations of apex meteors, except for the eccentricity distribution, which is 
more skewed toward e ~ 1 in our model. Finally, we propose an explanation for 
the long-standing problem in meteor science related to the relative strength of 
apex and helion/antihelion sources. As we show in detail, the observed trend, 
with the apex meteors being more prominent in observations of highly sensitive 
radars, can be related to orbital dynamics of particles released on the long-period 
orbits. 
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1. Introduction 

The Oort cloud is a roughly spherical cloud of comets (Oort 1950), which surrounds 
the solar system and extends to heliocentric distances larger than 100,000 AU. The Oort 
cloud is currently feeding comets into the inner solar system at a rate of about 12 comets 
with q < 3 AU per year with an active comet absolute magnitude H w < 10.9 (Wiegert & 
Tremaine 1999; q is the perihelion distance and H io is a distance-independent measure of the 
active comet brightness that includes the coma). The Oort Cloud Comets (hereafter OCCs) 
can be divided into two populations based on their dynamical histories: (1) dynamically new 
OCCs, which are on their first passage through the inner solar system and typically have 
a > 10,000 AU; and (2) returning OCCs which have previously passed through the inner 
solar system and typically have a < 10, 000 AU. 1 

The dynamical models of the orbital evolution of new OCCs into returning OCCs predict 
many times more returning comets than are observed (Wiegert & Tremaine 1999). This is 
the so-called “fading problem”, which is thought to be related to the physical evolution 
of OCCs. To resolve this problem, Levison et al. (2002, hereafter L02) proposed that 
OCCs must physically disrupt as they evolve inward from the Oort cloud. Specifically, 
L02 estimated that, when a OCC becomes inactive, it has only ~1% chance of becoming 
dormant, and ~99% chance of being disrupted. If more OCCs would become dormant, 
L02 argued, the modern surveys of near-Earth objects would discover a greater number of 
dormant OCCs passing through perihelion each year than they do. Strong thermal gradients, 
phase transitions, and volatile pressure buildup experienced by OCCs during their approaches 
to the Sun are thought to be responsible for disruptions. 

If these results are correct, the disrupted OCCs must be a prodigious source of dust 
particles and larger fragments that may further disintegrate. In particular, the dust pro- 
duction rate from OCC disruptions should be vastly larger than that of active OCCs. On 
the other hand, the smallest dust particles produced in these disruption events may be lost 
from the Solar System due to the effects of radiation pressure, while the large fragments 
should be dispersed over enormous radial distances. It is therefore not clear whether the 


1 See Dybczynski (2001) for an alternative definition. 
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disrupted OCCs can supply a significant amount of material to the inner zodiacal cloud, and 
whether they could represent a significant source of interplanetary dust particles accreted 
by the Earth. 

Given their large speeds relative to the Earth, the OCC particles plunging into the 
upper atmosphere could produce meteor phenomena and be detected by optical and radar 
meteor surveys (e.g., see Steel (1996) for a review). Here we consider the radar meteors. The 
modern meteor radar surveys produced vast datasets including millions of high quality orbits 
covering both the northern and southern hemispheres (e.g., Jones & Brown 1993, Taylor & 
Elford 1998, Galligan & Baggaley 2004, 2005, Janches et al. 2003, Janches & Chau 2005, 
Chau et al. 2007, Campbell-Brown 2008). Moreover, the sensitivity of meteor observations 
broke new grounds with routine measurements of meteor echoes using HPLA 2 radars such 
as the Arecibo (AO) radar, which is capable of detecting ~50-/rm particles down to ~20 km 
s^ 1 (e.g., Janches et al. 2003, 2006, Fentzke et al. 2009). It is natural to ask whether some 
of these observations can be linked to the particle populations from disrupted OCCs. 

Meteors are produced by small interplanetary particles, also known as the meteoroids, 
that interact with air molecules upon atmospheric entry. Based on meteor data, the me- 
teoroids can be divided into two groups: sporadic meteoroids and meteoroid streams. The 
meteoroid streams are prominent concentrations of particles with similar orbits (Whipple & 
Gossner 1949, Whipple 1951). They are thought to be produced by particles released by 
active and recently (< few thousand years ago) disrupted comets (e.g., Jenniskens 2008). 
Sporadic meteoroids are those particles that have evolved significantly from their parent 
body so that they are no longer easily linked to that parent, or to other meteoroids from the 
same parent. Notably, the time-integrated flux of meteors at Earth is dominated by about 
a factor of ~10 by sporadics (Jones & Brown 1993). 

The radiant distribution of sporadic meteors shows several concentrations on the sky 
known as the helion/antihelion, north/south apex, and north/south toroidal sources (e.g., 
Younger et al. 2009, and the references therein). Wiegert et al. (2009; hereafter W09) 
developed a dynamical model to explain these concentrations. Their main results concern 
the prominent helion/antihelion sources for which the particles released by Jupiter-Family 


2 High Power and Large Aperture. 



4 


Comets (JFCs) such as 2P/Encke provide the best match, in agreement with previous studies 
(e.g., Jones et al. 2001). As for the north/south apex source, W09 pointed out the potential 
importance of retrograde Halley-type comets (HTCs) such as 55P/Tempel-Tuttle (or an or- 
bitally similar lost comet). The case for the retrograde HTC particles is compelling, because 
three prominent retrograde HTCs, namely lP/Hallcy, 55P/Tempcl- Tuttle and 109P/Swift- 
Tuttle, all have associated streams, known as Orionids /rj Aquarids, Leonids and Perseids, 
respectively. The sporadic meteoroids with the north/south apex radiants can thus plausibly 
be a dynamically old component of HTC particles. 

Here we consider the possibility that at least some part of the meteoroid complex is 
produced by disrupting OCCs (L02, see also Jones et al. 2001). We study the effects 
of radiation pressure on particles released from the highly-eccentric OCC orbits, and their 
dynamical evolution under gravitational perturbations from planets and Poynting-Robertson 
(P-R) drag (see Section 2 for our model). We show that a significant contribution of OCC 
particles to the inner zodiacal cloud and meteor record is somewhat problematic, because 
most small OCC particles are blown out from the solar system by radiation pressure, while 
most large ones get scattered by planets and never make it into the inner solar system (Section 
3). Still, we find that there is a sweet spot at particle sizes ~100-300 /im. Our modeling work 
shows that the orbits and impact speeds of these intermediate-size OCC particles can match 
those derived from the meteor radar data for apex meteoroids. Furthermore, we find that 
the preponderance of fast apex meteors in HPLA radar observations (e.g., AO, ALTAIR, 
Jicamarca) can be linked to the competing effects of P-R drag and Jupiter perturbations, 
which act as a size filter on populations of the long-period meteoroids (Section 4). 


2. Model 

We studied the following sequence of events: (1) particles of different sizes were released 
from OCCs (Section 2.1), (2) their orbits evolved under the influence of gravitational and 
radiation forces (Section 2.2), (3) some particles were thermally or collisionally destroyed 
(Section 2.3), and (4) a small fraction of the initial particle population was accreted by 
Earth, producing meteors (Section 2.4). We describe our model for (l)-(4) below. 
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2.1. Initial Orbits 

According to Francis (2005), OCCs have d N(q) oc (1 + y/q) dq for q < 2 AU. For 
q > 2 AU, Francis’ result predicts d N(q) being flat or declining while we would expect the 
perihelion distribution to increase with q. It probably just shows that the distribution is 
not well constrained for q > 2 AU. We used d N(q) oc 2.41(g/2) 7 dq for q > 2 AU, with 
0 < 7 < 1. The initial values of q in our numerical integrations were set to be uniformly 
random between 0 and 5 AU, because particles starting with q > 5 AU do not reach 1 AU 
(see Section 3.1), where they could contribute to the Earth impact record. The results for 
different dN ( q ) were obtained by assigning the appropriate weight to particles starting with 
different g’s. 

Upon its release from a larger object a small particle will feel the effects of radiation 
pressure. These effects can be best described by replacing the mass of the Sun, m e , by 
m Q (l — fi), with P given by 

P = 5.7 xlO" 5 ^, (1) 

ps 

where radius s and density p of the particle are in cgs units. Pressure coefficient Q pr can be 
determined using the Mie theory (Burns et al. 1979). We set Q pr = 1, which corresponds to 
the geometrical optics limit, where s is much larger than the incident-light wavelength. We 
used particles with D = 2s = 10, 30, 100, 300, 1000 pm, which should cover the interesting 
range of sizes, and p = 2.0 g cm -3 . 

For large eccentricity e of the parent object and/or for large P, the released particle may 
become unbound and escape to interstellar space. To stay bound, the heliocentric distance, 
R , of the released particle must fulfill the following condition (e.g., Kresak 1976, Liou et al. 
1999): 

R > R* = 2pa . (2) 

This condition shows that all particles with P released at the orbit’s perihelion will be 
removed, if 2/3 > 1 — e. The new OCCs have 1 — e < 10~ 4 . It follows that particles produced 
by a new OCC near its perihelion will become unbound for sizes up to D ~ 1 cm. The usual 
near-perihelion activity of OCCs therefore cannot be a major source of small dust particles 
in the inner solar system. 
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Motivated by the L02 results, we now consider OCCs disruptions. Interestingly, obser- 
vations of the disruption events of comets show that there does not seem to be any correlation 
between the time of disruption and the orbital phase of the parent object. Many comets 
were seen to disrupt (or suffer outburst/splitting events) at large heliocentric distances. For 
example, 174P/Echeclus showed an outburst with R « 13 AU, more than 6 AU beyond its 
perihelion distance (Choi et al. 2006). It may therefore be possible that OCCs could disrupt 
at relatively large R and produce particles that, according to Eq. (2), will stay on bound 
orbits. 

We release particles with R > R* in our model. For example, a D — 1 00-// 1 i i particle 
with p = 2.0 g cm” 3 ejected from the parent comet with a = 10 3 AU will have (3 — 0.006 and 
R* = 12 AU. We thus release these particles with R > 12 AU. In addition, we only study 
particles ejected from orbits similar to those of the returning OCCs with a rs_/ 10 3 AU. We do 
not consider orbits with a > 10 4 AU, corresponding to the Oort spike 3 , because we believe 
it unlikely that disruptions could happen at the very large heliocentric distance implied by 
Eq. (2) for a > 10 4 AU. For example, a D — 300-/nn particle with p = 2.0 g cm” 3 , released 
from a parent orbit with a = 10 4 AU, would become unbound, unless R > 36 AU. 

The particle populations studied here have bound initial orbits. They represent only 
a fraction of all particles released from OCCs. This fraction, denoted /o in Section 4, is 
difficult to estimate, because we do not have a detailed understanding of the processes, 
and their dependence on R, that govern comet disruptions. We will return to this issue 
in Section 4. The initial distribution of orbital inclination vectors was set to be isotropic. 
To simplify things, we neglected the ejection velocities of dust particles from their parent 
objects (see Jones et al. 2001) and assumed that they will initially follow the parent comet’s 
orbit modihed by radiation pressure. 


3 The semimajor axis values of most OCCs are 10 4 <a<5x 10 4 AU, which is known as the Oort spike 

(e.g., Wiegert & Tremaine 1999). Comets in the spike are mostly dynamically new comets, on their first 
passage into the inner planetary system from the Oort cloud. A comet that passes through the planetary 
system receives a gravitational kick from the planets. The typical energy kick, Aa;, depends strongly on the 
perihelion distance of the comet’s orbit. According to Wiegert & Tremaine (1999), Aa: ~ 10” 3 AU” 1 for 
q < 6 AU, while comets in the Oort spike have x = 1/a < 10” 4 AU” 1 . Depending on the sign of the kick, 
they will either leave the planetary system on unbound orbit, never to return, or be thrown onto a more 
tightly bound orbit with a < 10 3 AU. 
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For reference, we also followed meteoroids from lP/Halley, 2P/Encke and 55P/Tempel- 
Tuttle. These comets were suggested to be important sources of the sporadic meteors by 
W09. The comet’s orbits were obtained from JPL fforizons site. Particles of different 
sizes were released from the parent orbits and tracked into future. We applied the same 
procedures/criteria to them that we used for the OCC particles. 


2.2. Orbit Integration 


The orbits of small particles in the interplanetary space are subject to gravitational 
perturbations of planets and radiation forces (Robertson 1937, Burns et al. 1979). The 
acceleration F due to radiation forces is 


F 


pG 


m 0 

R 2 




R 

R 


V 

c 


( 3 ) 


where R is the heliocentric position vector of particle, V is its velocity, G is the gravitational 
constant, m® is the mass of the Sun, c is the speed of light, and R = dR/dt. The acceleration 
(3) consists of the radiation pressure and the velocity-dependent P-R term. Parameter P is 
related to the radiation pressure coefficient Q pi by Eq. (1). 

The particle orbits were numerically integrated with the swift_rmvs3 code (Levison & 
Duncan 1994), which is an efficient implementation of the Wisdom- Holman map (Wisdom & 
Holman 1991) and which, in addition, can deal with close encounters between particles and 
planets. The radiation pressure and drag forces were inserted into the Keplerian and kick 
parts of the integrator, respectively. The change to the Keplerian part was trivially done by 
substituting m® by m®(l — P). The swift_rmvs3 integrator is stable even for near-parabolic 
orbits, and thus well suited for the integrations that we carried out here. 

The code tracks the orbital evolution of a particle that revolves around the Sun and is 
subject to the gravitational perturbations of seven planets (Venus to Neptune; the mass of 
Mercury was added to the Sun) until the particle impacts a planet, is ejected from the Solar 
System or evolves to within 0.05 AU from the Sun. We removed particles that evolved to 
R < 0.05 AU, because the orbital period for R < 0.05 AU was not properly resolved by our 
1-day integration timestep. 



Several thousand particles were followed for each D. Their orbital elements were defined 
with respect to the barycenter of the solar system. The barycentric elements are similar to 
the heliocentric elements for R < 5 AU, but differ for large R, where the Sun’s orbital speed 
about the barycenter of the solar system is not negligible relative to the orbital speed of a 
particle. 

All orbits were followed from the present epoch into the future. Each particle’s orbital 
elements were stored at 10 3 yr intervals. We used the output to construct a steady state 
distribution of OCC particles in the inner solar system. This approach differs from that of 
Nesvorny et al. (2006) and W09, who started particles at many different past epochs and used 
these integrations to determine the present distribution of particles. The two distributions 
are expected to be slightly different expressing mainly the difference between the present 
configuration of planets, with each planet having a specific secular phase, and the time- 
averaged system, where all phases are mixed. Since this difference is small, however, we can 
use the steady state distribution, which is easier to obtain, as a reasonable approximation. 


2.3. Physical Effects 

Solar system micrometeoroids can be destroyed by collisions with other particles and by 
solar heating that can lead to sublimation and vaporisation of minerals. Here we describe 
how we parametrize these processes in our model. 


2.3.1. Thermal Destruction 

Thermal alteration of grains in the interplanetary grains is a complex process. The OCC 
particles evolving into the inner solar system will first loose their volatile ices, which will 
rapidly sublimate once the grains are heated to a critical temperature. We do not model the 
volatile loss in detail. Instead, we crudely assume that the grains have lost ~50% of their 
mass/volume when reaching R < 5 AU. We do not include the orbital effects of mass loss 
in orbital modeling, because it should produce only a relatively small perturbation on orbits 
for large particles that we consider here. The remaining grains will be primarily composed 
from amorphous silicates and will survive down to very small R. 
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According to Duschl et al. (1996), silicates are thermally altered at temperatures T ~ 
900-1600 K, and start to vaporise for T > 1600 K. As an example of thermal alteration, 
Kasuga et al. (2006, see also Capek & Borovicka 2009) studied the thermal desorption 
of alkali minerals and concluded that micrometeoroids should show evidence of thermal 
desorption of metals, Na in particular, for q < 0.1 AU. Following Moro-Martfn & Malhotra 
(2002), Kessler-Silacci et al. (2007) and others, we adopt a simple criterion for the silicate 
grain destruction. We will assume that they are destroyed when the grain temperature 
reaches T = 1500 K. 

The temperature of a small, fast spinning grain in interplanetary space is set by an 
equilibrium between the absorbed and re-radiated energy fluxes. While the absorbed flux is 
a simple function of the particle’s size, albedo and heliocentric distance, the re-radiated flux 
depends on the particle’s emissivity, which in turn is a function of particle’s size, shape and 
material properties. Using the optical constants of amorphous pyroxene of approximately 
cosmic composition (Henning & Mutschke 1997), we find that a dark D > 100-/mi grain 
at R has the equilibrium temperature within 10 K of a black body, T(R ) ~ 275 / \fR K. 
According to our simple destruction criterion, T(R) > 1500 K, the silicate grains should 
thus be removed when reaching R < 0.03 AU. On the other hand, the smallest particles 
considered in this work, D — 10 // m , will reach T(R) = 1500 K for R ~ 0.05 AU. Thus, we 
opted for using a very simple (and conservative) criterion where particles of all sizes were 
destroyed, and not considered for statistics, if they ever reached R < 0.05 AU. Note that, by 
design, this limit is the same as the one imposed by the integration timestep (Section 2.2). 


2.3.2. Disruptive Collisions 

The collisional lifetime of meteoroids, r co ii, was taken from Grim et al. (1985; hereafter 
G85). It was assumed to be a function of particle mass, m, and orbital parameters, mainly 
a and e. We neglected the effect of orbital inclination on r co u, because the results discussed 
in Steel & Elford (1986) suggest that the inclinations should affect r co u only up to a factor of 
~2-5, which is not overly significant in the context of our work. We assumed that the mass 
and orbital dependencies of r co ii can be decoupled, so that 


T co n (m, a, e) = $ (m) T (a, e) , 


( 4 ) 
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where $(m) and T(a, e) are discussed below. 

As for <h(m), we used the G85’s model based on measurements of various spacecraft and 
Earth-based detectors. We found that the model can be approximated by a simple empirical 
fit. Specifically, between 10~ 64 g and 10 1-2 g, we adopted the quadratic function 

log (m) = c 2 (logm) 2 + (logm) + c 0 , (5) 


with (c 0 ,ci,c 2 ) = (4.021,0.300,0.083), where the values of c 0 , C\ and c 2 were set to fit 
the G85’s collision lifetime for circular orbits at the reference distance i? 0 — 1 AU. A linear 
relation between log (m) and log m was used outside the quoted mass range to approximate 
the G85’s model down to m = 10 -8 g. $(m) has a minimum for m ~ 0.01 g, corresponding 
to s ~ 1 mm for p = 2 g cm -3 (Fig. 1). The collisional lifetime of ~1 mm particles in the 
G85’s model is very short, roughly 5,000 yr at 1 AU. 


T(a, e) is assumed to drop as a power-law with R. From Eq. (18) in G85 we have 


( R\ a v{R) 

\Ro J ^Vii'c ( R) 


( 6 ) 


where v(R ) and v C i rc (R) are the particle and circular speeds at R , respectively, and a ~ 1.8. 
The velocity-dependent factor provides an appropriate scaling of r co ii for eccentric orbits. 


Averaging Eq. (6) over a Keplerian orbit with semimajor axis a and eccentricity e, we 
obtain 


T (a, e) = 


Rr 


J(e) • 


where 


J(e) = 


du 


(1 


e cos u 


\l+a 


7T 


( 7 ) 


( 8 ) 


a/1 + ecos u 

can be written as a series in e 2 with good convergence. Note that J ~ 1 for e ~ 0, as 
required in Eq. (7), but can become 3>1 for very eccentric orbits. 


The G85’s model was calibrated to match the impact fluxes of particles as measured 
prior to 1985. The more recent measurements indicate lower fluxes (e.g., Dikarev et al. 
2005, Drolshagen et al. 2008). Also, to estimate r co n, assumptions needed to be made in 
G85 about the strength of particles, and their impact speeds. As a result, r co n proposed by 
G85 may have a significant uncertainty. To test different possibilities, we introduced two free 
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parameters in our model, Si and S 2 , that were used to shift the <3>(m) function in logm and 
log<f>, respectively (as indicated by arrows in Fig. 1). For example, the positive S 2 values 
increase t co u relative to the standard G85’s model, as expected, for example, if particles were 
stronger than assumed in G85, or if the fluxes were lower. 

Collisional disruption of particles was taken into account during processing the output 
from the numerical integration described in Section 2.2. To account for the stochastic nature 
of breakups, we determined the break-up probability p co \\ = 1 — exp(— h/r co n), where h = 
1000 yr is the output interval, and r co u was computed individually for each particle’s orbit. 
The code then generated a random number 0 < x < 1, and eliminated the particle if x < p co \\. 

We caution that our procedure does not take into account the small debris fragments 
that are generated by disruptions of larger particles. Instead, all fragments are removed from 
the system. This is an important approximation, whose validity needs to be tested in the 
future. 


2.4. Model for Meteor Radar Observations 

We used the Opik theory (Opik 1951) to estimate the expected terrestrial accretion rate 
of OCC particles in our model. Wetherill (1967), and later Greenberg (1982), improved the 
theory by extending it more rigorously to the case of two eccentric orbits. Here we used the 
Fortran program written by W. F. Bottke (see, e.g., Bottke et al. 1994), which employs the 
Greenberg’s formalism. 

We modified the code to compute the radiants of the impacting bodies. The radiants 
were expressed in the coordinate system, where longitude l was measured from the Earth’s 
apex in counter-clockwise direction along the Earth’s orbit, and latitude b was measured 
relative to the Earth’s orbital plane. Note that our definition of longitude is different from 
the one more commonly used for radar meteors, where the longitude is measured from the 
helion direction. The radiants were calculated before the effects of gravitational focusing 
were applied. 

The longitude and latitude values of radiants were binned into 1 deg 2 area segments. 
For each radiant bin, the code gives information about the distribution of geocentric impact 
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speeds, v g , and heliocentric orbits prior to the impact, as defined by a, e and i. Here, v g 
is defined as v g = (v^ + v g SC ) 1//2 , where is the relative velocity ‘in infinity’ and v esc is 
the escape speed from the Earth’s surface. The orbital elements, on the other hand, are the 
orbital elements that the particle would have in absence of the gravitational focusing by the 
Earth. The radiant distribution, v g , a, e and i will be compared to meteor radar observations 
in Section 3.3. 


To compare our model with observations, we need to include the meteor radar detection 
efficiency. This is a difficult problem because the meteor phenomenon itself and radar detec- 
tion of it involve complex physics. For example, the specular meteor radars (SMRs), such as 
the Canadian Meteor Orbit Radar (CMOR; Campbell-Brown 2008) and Advanced Meteor 
Orbit Radar (AMOR; Galligan & Baggaley 2004, 2005), detect the specular reflection of the 
meteor trail (the plasma formed by the meteoroid’s passage). The meteoroid velocities are 
then derived from the detection of the Fresnel diffraction patterns of the developing trail, or 
are determined by measuring the time of flight between stations. 


The detection efficiency of a meteor should mainly be a function of the particle size 
and speed, but it also depends on a number of other parameters discussed, for example, 
in Janches et al. (2008). Following W09, we opt for a simple parametrization of radar 
sensitivity function, where the detection is represented by an ionization function 

/ \ 3.5 


I(m,v g ) 


m I Vg 

10~ 4 g \30 krri/s 


( 9 ) 


All meteors with I(m,v g ) > I* are assumed to be detected in our model, while all meteors 
with I(m,v g ) < I* are not detected (see Fentzke et al. (2009) for a similar method applied 
to head echo radars). The ionization cutoff I* is different for different SMRs. For example, 
I* ~ 1 for CMOR (Campbell-Brown 2008) and I* ~ 0.001-0.01 for AMOR (Galligan & 
Baggaley 2004, 2005). For reference, an OCC particle with v g = 60 km s -1 and m = 10 -5 g, 
corresponding to s ~ 100 /jm, will have I(m,v g ) ~ 0.1, i.e., a value intermediate between 
the two thresholds. These meteoroids would thus be detected by AMOR, according to our 
definition, but not by CMOR. We will discuss these issues in more detail in Section 3.3. 
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3. Results 

3.1. Orbital Evolution of OCC Particles 

Wyatt & Whipple (1950) identified the following constant of motion of P-R drag 

c = < 10 > 

where a and e are the particle’s semimajor axis and eccentricity (see also Breiter & Jackson 
1998). This constant is independent of the particle properties such as its size. The orbit 
path of any particle in a and e can thus be obtained by calculating C for the initial orbit, 
and requiring that Eq. (10) holds at all times (Fig. 2a). For e ~ 1, the orbit trajectories 
follow the lines of constant q, because C ~ 2q for e ~ 1. The semimajor axis of orbits shrinks 
till reaching a value only several times larger than q. At that point, a more familiar form of 
P-R drag takes place with both a and q converging to zero. 

The timescale of orbital evolution is as follows. For an initial orbit with a, e and 
q = a(l — e), the total time of fall to the Sun is (Wyatt & Whipple 1950) 

r M , = 20 Myr (^) (j^s) (uu)' f(e) ' (U) 

where 

( 1 4- e ) 2 r e n 3 / 5 

/« = 1 d "(Tvy)W ' < 12 > 

Figure 2b shows Tf a n for orbits relevant to OCC particles. 

Interestingly, the timescale for large a can be relatively short if q is small. For example, 
an OCC particle with D = 100 /jrn, a = 10 4 AU and q — 1 AU takes about 20 My to fall to 
the Sun. This is only about 20 orbital periods for a = 10 4 AU. According to Eq. (11), Tf a n 
scales linearly with particle size. Thus, a D — 10-/jm particle with the same initial orbit has 
Thai = 2 My. 

One important aspect that cannot be captured by the analytical results discussed above 
is the effect of planetary perturbations on drifting orbits of OCC particles. To evaluate this 
effect, as described in Section 2.2, we numerically integrated the orbits of OCC particles as 
they evolve from large a and interact with the planets. 
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Figure 3 shows the orbital history of a particle, whose orbit evolved all the way down 
into the inner solar system. Initially, the particle’s semimajor axis underwent a random 
walk caused by indirect planetary perturbations, mainly from Jupiter, during each perihelion 
passage (Wiegert & Tremaine 1999). Then, at time t ~ 5 x 10 5 yr, a single perihelion passage 
produced a significant drop of a from 150 to 40 AU, where the particle started to interact 
with the exterior mean motion resonances with Neptune. The following evolution was mainly 
controlled by P-R drag. Eventually, the orbit decoupled from Jupiter, reached a ~ 1 AU, 
and kept shrinking further toward the Sun, where the particle was thermally destroyed. The 
large oscillations of q for 0.9 < t < 1.7 Myr, correlated with those in i, were produced by 
Kozai dynamics (e.g., Kozai 1962). 

The orbit history shown in Fig. 3 is typical for an OCC particle that is able to make it 
into the inner solar system. These particles, however, represent a relatively small fraction of 
the initial population, with most particles being ejected from the solar system by planetary 
perturbations. Using d N(q) with 0 < 7 < 1 (Section 2.1), a ~ 10 3 AU and isotropic 
distribution of inclination vectors, roughly 0.8-1. 5 % of particles with D = 100 /jui evolve 
down into the inner solar system and decouple from Jupiter (as defined by the aphelion 
distance Q = o(l + e) < 4 AU), without being disrupted by an impact (with standard 
Si = S 2 = 0; Section 2.3) or ever having q < 0.05 AU (to avoid sublimation). These 
particles can potentially be important for the terrestrial impact record and interpretation 
of the meteor radar data. The bulk of OCC particles that do not reach a ~ 1 AU, do not 
significantly contribute to the Earth’s accretion, because these particles spend most of their 
lifetimes at R ^ 1 AU. 

The fraction of OCC particles reaching a ~ 1 AU, fi, normalized to the number of 
particles whose orbits were initially bound (see Section 2.1), is sensitive to particle size. For 
D = 10 /jui, d N(q) described in Section 2.1, and initial a ~ 10 3 AU, fi ~ 0.15 — 0.2. For 
D = 300 /irn, on the other hand, /i~2x 10~ 3 . Moreover, for D = 1 mm, only one particle 
out of the total of 5000 reached a ~ 1 AU. This trend, with the larger particles having 
progressively smaller /j values, has interesting implications for the observations of sporadic 
meteors (Section 4). 

The above estimates used r co u as described in Section 2.3.2. Collisional disruption, 
however, turned out to have only a modest effect for the standard G85’s r co u (Ai = S 2 = 0) 



and the particle sizes considered here. For example, only ~2% of particles with D = 100 /mi 
that reached a ~ 1 AU in our numerical integration have disrupted prior to decoupling from 
Jupiter, as detected in post-processing of the integration output, with the standard G85’s 
T co u. A great majority of particles with D = 300 /mi also survived. The effect of disruptive 
collisions becomes more significant for D ~ 1 mm, for which r co n is significantly shorter 
(Section 2.3.2) and P-R drag is weaker. 

The thermal effects discussed in Section 2.3.1 turned out to be very important for all 
particle sizes considered here. For example, 75 out of 122 particles (i.e., over 60 %) with 
D = 100 /mi that ever reached Q < 4 AU, previously had q < 0.05 AU, which is our crude 
threshold for the thermal destruction of particles. Also, 68 out of 74 D = 300 /mi particles 
(over 90 %) reaching Q < 4 AU previously had q < 0.05 AU. The thermally destroyed 
particles are removed and the orbital histories after their disruption are not used for our 
analysis. 

These fractions are a direct consequence of the relative importance of planetary pertur- 
bations and P-R drag on particles with different D and q. The orbital evolution of a large 
particle on the OCC-like orbit is primarily controlled by planetary perturbations. Sooner 
or later, the planets will eject the particle from the solar system, unless the orbit shrinks 
and decouples from Jupiter. To achieve this, q of the particle’s orbit must be very low, 
so that the P-R drag timescale is short (see Fig. 2b). But if q is low, it may easily drop 
below q < 0.05 AU, where the particle is removed, thus explaining why most large particles 
reaching Q < 4 AU previously have q < 0.05 AU. 

Interestingly, the fraction of particles with D ~ 100-300 /mi reaching a ~ 1 AU is not 
overly sensitive to the initial perihelion distance, as far as q < 5 AU. This is because it is 
more likely to decouple if q is low, because P-R drag is stronger, but this trend is nearly 
canceled, because particles with very low q tend to drop below q = 0.05 AU and sublimate 
before they can decouple. Particles with D ~ 100-300 /mi and q > 5 AU, on the other hand, 
tend to have very long P-R drag timescales (e.g., Tf a n > 50 Myr for initial a = 10 3 AU; Fig. 
2b), and are scattered by planets from the solar system. 



3.2. Orbital and Spatial Distributions 


Here we discuss the expected distribution of OCC particles in the inner solar system. 
Figure 4 shows the number density of OCC particles as a function of R. As expected, the 
particle density increases toward the Sun. The radial distribution of D = 100 /in i particles 
can be approximated by a power law, d N(R) oc R~ a dR, where a ~ 1.5 for R < 5 AU, 
and a ~ 2.0 for R > 10 AU. Both these radial dependencies are significantly steeper than 
d N(R) oc R~ 1 dR, expected for distribution of particles on nearly-circular orbits (see, e.g., 
Dermott et al. 2001). 

The relatively steep radial distribution is a consequence of Keplerian motion of particles 
with e ~ 1. For R < a, the time spent by a particle on the Keplerian orbit between R and 
R + dR is dt oc R 0 5 dR. This leads to d N(R) oc R~ 15 dR, if the expected number of particles, 
which is proportional to dt, is divided by volume 47ri? 2 di?. For R > a, on the other hand, 
appropriate for large radial distances, dt ~ const, and d N(R) oc R~ 2 dR. Thus, a ~ 1.5-2 is 
expected for e ~ 1 (see also Liou et al. 1999). 

We find that the number density of OCC particles at R ~ 1 AU is mainly contributed 
by the particles that orbitally decoupled from Jupiter. This shows the importance of orbital 
decoupling for the distribution of OCC particles in the inner solar system, and their accretion 
by the Earth. Specifically, most OCC particles accreted by the Earth are expected to have 
a(l + e) < 4 AU despite the fact that their orbits started with a > 10 3 AU. 

Figure 5 shows the distributions of orbital elements for OCC particles with D = 100 /jrn 
and R < 5 AU. The semimajor axis distribution, dA^(a), has a broad maximum centered at 
0.5-2 AU with a ~ 2 AU being the most common. dN(a) decreases toward larger a, because 
particles with a > 3 AU are coupled to Jupiter, have short dynamical lifetimes, and do not 
spend much time at R < 5 AU. The peak of dN(a) at a = 6-7 AU is produced by orbits in 
the exterior mean motion resonances with Jupiter, which prolong the dynamical lifetime of 
particles by phase protecting them against encounters with Jupiter (Liou et al. 1999). 

The eccentricity distribution dN(e) increases toward e ~ 1, which is expected because 
all orbits started with e > 0.995. The tail extending to the moderate and low eccentricity 
values is due to the dynamically long-lived particles, whose orbits decouple from Jupiter and 
become circularized by P-R drag. 
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The inclination distribution is also interesting as it significantly deviates from the initial 
distribution with d N(i) oc sin(i)di (Fig. 5c). The retrograde orbits are more common than 
the prograde ones. The preference for retrograde orbits is probably caused by gravitational 
perturbations from Jupiter that are more effective on prograde orbits, because the encounter 
speeds are lower, and thus Ah’s are larger. The prograde particles should therefore have 
shorter dynamical lifetimes than the retrograde particles, which would explain their relative 
paucity in a steady state dN(i) for R < 5 AU. Interestingly, however, the particles that 
decouple from Jupiter and reach a ~ 1 AU by P-R drag do not show a strong preference for 
retrograde orbits. 

In addition, the steady state dN(i) lacks orbits with i ~ 90°. We believe that this 
is a consequence of Kozai dynamics (e.g., Kozai 1962). It is well known that the initially 
near-polar orbits will suffer large oscillations of q driven by variations of the orbital angular 
momentum vector. Most of these orbits can therefore reach very low q values, where the 
particles will be destroyed by thermal effects. Indeed, by studying the orbital histories of 
particles that started with i ~ 90°, we found that most of these orbits had q < 0.05 AU 
prior to reaching Q < 4 AU. 

The steady state distribution of orbital elements of D = 300-//m particles is similar to 
the one discussed above. Instead of having a single peak at a = 6-7 AU, however, dN(a) 
for a > 5 AU is more irregular showing many peaks and dips. Apparently, since the larger 
D = 300-yum particles drift more slowly by P-R drag, they are more susceptible to capture 
in a large number of resonances (Liou et al. 1999). The second difference concerns dA(e), 
which is slightly more clumped toward e ~ 1 for D = 300 //in than for D = 100 //in. 


3.3. Radiants and Orbits of Particles Accreted by the Earth 

The radiants of OCC particles are located near the Earth’s apex (Fig. 6). This is 
logical because the retrograde OCC particles, which come from the apex direction, have 
much higher velocity relative to the Earth (~60 km s -1 ) than particles on prograde orbits. 
The retrograde particles therefore also have, according to the usual nav rule, rather large 
impact probabilities with the Earth. In addition, the ionization cutoff used here (Eq. 9) 
poses a rather strict limit on the mass of meteoroids that can be detected by SMRs at low 
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speeds. For example, a D = 100-/mi particle with m = 10 -6 g and v g = 30 km s~ x has 
/ = 0.01, which is near or slightly above the detection limit of AMOR, and way below the 
detection limit of CMOR. 

The model radiants form the south and north apex concentrations, just as observed 
(e.g., Jones & Brown 1993, Galligan & Baggaley 2005, Campbell-Brown 2008). The lack 
of radiants within ~10° about the ecliptic is due to near absence of OCC meteoroids with 
i ~ 180° (see Fig. 5c). The lack of radiants with b > 50° (or b < —50°) is the consequence 
of the inclination distribution shown in Fig. 5c that is depleted in orbits with i ~ 90°. 

Most north (south) radiants fall into an area on the north (south) hemisphere that has 
the characteristic triangular or half-disk shape. For D = 100 //rn, the centers of radiant 
concentrations are at b ~ ±20°. For D = 300 //rn, the centers are at b ~ ±25°. This 
reflects the differences in d N{i) between the populations of particles with D = 100 /jrn and 
D = 300 /jrn that we obtained in the model. For comparison, observations indicate that 
b ~ ±15° (e.g., Galligan & Baggaley 2005, Chau et al. 2007, Campbell- Brown 2008). In 
addition, the apex sources that we obtained in our model tend to be more stretched in both 
l and b than the observed ones. 

While more modeling work will be needed to test things with a better statistic, the 
issues discussed above may indicate that a better fit to observations could be obtained if 
the retrograde source had inclinations closer to 180° than the bulk of OCC particles with 
100° < i < 160°. It is not clear how this could be achieved by tweaking the parameters 
of our model. Instead, clues such as these seem to highlight the importance of known 
HTCs. Indeed, the two prominent active HTCs, lP/Halley and 55P/Tempel-Tuttlc, both 
have i ~ 162°. They would therefore be expected to produce apex sources closer to the 
ecliptic than the bulk of retrograde OCC particles (Fig. 7). 

An interesting feature in Fig. 6b is the presence of a ring that stretches to ±60° in 
longitude and latitude. A similar ring has been noted in Campbell-Brown (2008), who 
suggested that the region inside the ring can be depleted in meteor radiants, except for apex 
sources, because the retrograde meteoroids with radiants inside the ring would have shorter 
collisional lifetimes. Our simple collisional model cannot reproduce this effect, because r co ii 
is independent of i (see Section 2.3.2). In addition, W09 suggested that the origin of the 
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ring can be traced back to Kozai dynamics, which confines the allowed radiants of particles 
on high-inclination orbits. Here we confirm the W09’s result by isolating particles that 
contribute to the ring, and checking on their orbital behavior. 

The impact speed of OCC particles peaks at v g ~ 60 km s _1 (Fig. 8 for D = 100 /mi and 
Fig. 9 for D = 300 /mi), which is a nice match to observations of apex meteors (Galligan & 
Baggaley 2005, Campbell- Brown 2008). 4 The width of the peak also looks good (cf. Fig. 13 
in Campbell-Brown 2008). In comparison, using a population of meteoroids from lP/Halley 
and 55P/Tempel-Tuttle, W09 obtained a peak at v g ~ 70 km s -1 , which is expected for large 
particles that have not evolved far from their parent comet orbit. 

Figures 8 and 9 also show the model distributions of orbital elements of OCC meteoroids. 
Distributions dlV(a), now heavily weighted by the collision probability (cf. Fig. 5), peak 
at a ~ 1 AU and show a tail extending to a > 2 AU. The meteoroids in the peak have 
orbits that have strongly evolved by P-R drag. The model distributions for D = 100 /mi 
and D = 300 /mi are similar, except for a few wiggles produced by statistical fluctuations. 
Both provide a good match to the observed semimajor axis distribution of apex meteors (cf. 
Fig. 11 in Galligan & Baggaley 2005). 

According to our model, most apex meteors should have inclinations between i ~ 100° 
and 180° (Figs. 8c and 9c), which is also the range indicated by observations. Unlike d N(i) 
measured by radars, which shows a broad peak centered at i ~ 150-160°, our model dN(i) 
is more spread and noisy. While part of this discrepancy could be blamed on insufficient 
statistics in our model, it may also point to a more fundamental problem. 

The eccentricity distribution is puzzling. For both D = 100 /mi and D = 300 /mi, 
we obtained dIV(e) that raises toward e ~ 1. This trend is slightly more pronounced for 
D = 300 /mi (Fig. 9d) than for D = 100 /mi (Fig. 8d). W09, using selected HTCs for parent 
bodies of apex meteors, obtained diV(e) that also peaked toward e ~ 1. In contrast, the 
observed apex meteors have nearly flat dIV(e) at 0.2 < e < 1, and show a slight depletion 
for e < 0.1 (Fig. 13 in Campbell-Brown 2008). 


4 Note that observations by High Power and Large Aperture (HPLA) radars such as AO or ALTAIR 
measure the apex peak speed at ~ 55 km s _1 . 
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The cause of these differences is unclear. To obtain lower values of e in our model, 
the orbits would need to become more circularized by P-R drag before arriving to a ~ 1 
AU. This could be achieved, for example, if more weight is given to particles starting with 
q > 1 AU. We confirm this by using 7 > 1 , but a detailed match to the observed eccentricity 
distribution remains elusive. A detailed analysis of this problem is left for future studies. 


4. Relative Importance of Helion/Antihelion and Apex Sources 

The observations of sporadic meteors show that the relative importance of hclion/anti- 
hclion and apex sources depends on the sensitivity of the radar that is used to carry out 
such observations. The less sensitive SMRs with I* ~ 1, such as the Harvard Radar Meteor 
Project (HRMP; Jones & Brown 1993, Taylor & Elford 1998) or CMOR, detect ~3-10 times 
more helion/antihelion meteors then apex meteors (see, e.g., Campbell-Brown (2009) for 
comparison of different radars). The more sensitive radars, such as AMOR with I* ~ 0.001- 
0.01, on the other hand, detect a relatively larger number of apex meteors. Finally, the apex 
meteors are predominant in observations by the highly-sensitive AO radar (e.g., Janches et 
al. 2003), because of their ability to detect small particles (Fentzke & Janches 2008, Fentzke 
et al. 2008). 

This trend can be explained if the size frequency distribution (SFD) of apex meteoroids 
is steeper (i.e. , if the number of meteoroids more sharply increases with decreasing size) 
than that of the helion/antihelion meteoroids, because radars that are capable of detecting 
smaller meteoroids would then be expected to see many more apex meteors (Fentzke & 
Janches 2008). For example, the initial SFD of particles produced by disrupted OCCs 
(or HTCs) could be steeper than the one produced by the sources of the helion/antihelion 
meteoroids (presumably active and disrupted JFCs; W09, Nesvorny et al. 2010). While this 
is a possibility that cannot be ruled out by the existing data, there are also no indications 
that this might be true (McDonnell et al. 1987, Grim et al. 2001, Green et al. 2004, Reach 
et al. 2007). 

We propose that the predominance of apex meteors in AO observations is caused by 
orbital dynamics of particles released from OCCs (or HTCs). Let us assume that the initial 
SFD of particles released from OCCs (or HTCs) is d N 0 (D) = N 0 D^^dD. The SFD of 



21 


meteoroids accreted by the Earth will then be 

d N(D) = foPiNoD^dD , (13) 

where / 0 is the fraction of particles that remain on bound orbits (see Section 2.1), and P x 
is the impact probability of these particles on the Earth. Factor f 0 expresses the removal 
of small particles by radiation pressure (see Section 2.1). Thus, the size dependence of /o is 
such that it cannot increase the number of small particles relative to the large ones. 5 

Table 1 lists P[ and Voo for various sources. For example, the OCC particles with 
D — 30 /jin and D = 300 /jrn have Pi = 2 x 10 6 and Pi = 5 x 10 7 , respectively. This 
indicates that dN(D) should have a steeper slope than dN 0 (D). Specifically, if dN 0 (D ) 
can be approximated by D~^dD for D ~ 30-300 pi n, where ( is a constant, we find that 
d N(D) oc D~^ +s MD, where 5 ~ 0.6 for OCCs. This estimate was obtained with 0 < 7 < 1, 
S\ ~ S 2 ~ 0, q* = 0.05 AU and initial a ~ 10 3 AU. The results for S 2 = 1, applicable if r co u 
were ~10 times longer than in G85, and for HTCs are similar. 

For comparison, W09 suggested that comet 2P/Encke (or an orbitally similar lost comet) 
is the main source of helion/antihclion meteors. We find that particles released from comet 
2P/Encke have Pi — 6 x 10 -5 for D = 30 pm and Pi = 2.3 x 10 -4 for D = 300 pm (Table 1). 
If these estimates are representative for the sources of the helion/antihclion meteors, they 
suggest that the slope of dN(D) should be shallower than that of cUVo(P) (5 ~ —0.6). 

We therefore find that dN(D) of apex meteoroids is expected to be steeper than that of 
helion/antihclion meteoroids, even if the initial SFD of particles released from the respective 
sources -JFCs and OCCs/HTCs- were similar. This effect is produced by orbital dynamics 
of particles starting on different initial orbits. As we discussed in Section 3.1, the Earth 
impact record of OCC particles is mainly contributed by those particle that decouple from 
Jupiter. Since the decoupling efficiency, described by factor f\ in Section 3.1, ramps up 
toward smaller D, the population of small OCC (or HTC) particles is enhanced, relative to 
large ones. This effect is weaker for JFC meteoroids, for which the correlation between 
and Pi is not as extreme. In the JFC case, the population of small particles accreted by the 


5 Factor fo{D) could presumably be approximated by a step function with fo{D) = 0 for D < D* and 

fo(D ) = 1 for D > D * , where D* is the critical diameter implied by Eq. (2). 



22 


Earth is suppressed by their short P-R drag timescale, and consequently, lower Pj. 

The magnitude of the SFD effects discussed above is just right to explain observations. 
If the ionization threshold I* of AMOR is ~100 times lower than that of CMOR/HRMP, 
these more sensitive instruments are expected to detect meteoroids that are ~5 times smaller 
in size. If they detect 3-10 times more apex meteors than the helion/antihelion meteors (e.g., 
Campbell-Brown 2008), this would suggest that the SFD slope index difference between apex 
and helion/antihelion meteors is 6 ~ 0.7-1. 4. We found 5 ~ 1.2 above, in a good agreement 
with observations. Figure 10 illustrates the relative strength of sporadic meteor sources 
expected from our model. 


5. Conclusions 

We found that only a very small fraction, f\ < 10 -4 , of D > 1-mm OCC particles can 
ever make it into the inner solar system. The relevance of these very large OCC particles 
to observations of sporadic meteors is therefore not obvious. The situation looks more 
favorable for OCC particles with D ~ 100-300 /jrn. These particles should survive the effects 
of radiation pressure, if released from returning OCCs at R > 4-12 AU. Moreover, about 
0.2- 1.5% avoid being collisionally disrupted or thermally destroyed, decouple from Jupiter, 
and finally spiral down to a ~ 1 AU, where their Earth-impact probability is increased by 
orders of magnitude. 

We estimated that the overall probability of Earth impact per one particle released on 
bound orbit from the returning OCC is Pj ~ 0.5-1 x 10 6 for D ~ 100-300 /jrn. This is 
50-80 times lower than P* expected for particles released from HTCs such as lP/Halley and 
55P/Tempel-Tuttle, and 200-400 times lower than P* expected for JFCs such as lP/Encke. 
The OCC particles will therefore significantly contribute to the sporadic meteor complex 
only if the mass of material produced by disrupting OCCs is large enough to compensate 
for these factors. From L02, we can roughly estimate that ~5 returning OCCs disrupt per 
year producing the mass input of perhaps as much as ~ 10 18 g yr _1 , or 3 x 10 8 kg s _1 . 
Only a small fraction of this mass will end in bound particles with D ~ 100-300 /jrn. For 
comparison, the active JFCs produce rx/ 300 kg s 1 (Reach et al. 2007). 
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We found that the SFD of apex meteoroids, presumably starting on highly-eccentric 
orbits, is expected to be steeper than those of helion/antihelion meteoroids, even if their 
initial SFDs were similar. The steepening of the SFD slope of apex meteoroids results from 
the efficiency with which OCC/HTC meteoroids of different sizes decouple from Jupiter. 
This result has interesting implications for observations of sporadic meteors, because it can 
explain why the north/south apex sources are more represented in observations of highly 
sensitive radars that are capable of detecting smaller meteoroids. 
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with us his Opik code, M. Campbell-Brown for useful discussions, and Tadeusz Jopek for a 
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km s -1 


OCCs 


10 

4.6 (4.6) 

52 (52) 

30 

1.9 (1.9) 

47 (47) 

100 

1.0 (1.4) 

46 (46) 

300 

0.5 (0.6) 

55 (54) 

1000 

0.2 (0.3) 

58 (58) 


lP/Halley (HTC) 

10 

57 (57) 

60 (60) 

30 

100 (110) 

59 (59) 

100 

80 (110) 

59 (59) 

300 

24 (46) 

65 (61) 

1000 

17 (22) 

67 (67) 

55P/Tempel- Tuttle (HTC) 

10 

53 (53) 

60 (60) 

30 

120 (120) 

59 (59) 

100 

80 (120) 

60 (59) 

300 

27 (51) 

66 (61) 

1000 

16 (23) 

68 (66) 


2P/Encke (JFC) 

10 

17 (17) 

18.0 (18.0) 

30 

60 (60) 

18.5 (18.5) 

100 

210 (220) 

18.7 (18.4) 

300 

230 (340) 

23.6 (22.0) 

1000 

120 (300) 

30.5 (27.0) 


Table 1: The Earth impact probability, Pi, and mean impact speed, {voo), of particles released 
on different orbits. For each particle’s diameter, D, we give our best estimate values for the 
standard G85’s t co u (Si — S 2 — 0), and for Si = 0 and S 2 = 1 (values in parenthesis). 
The longer collisional lifetime in the later case leads to the larger P t values. The effect of 
disruptive collisions is significant for D > 300 // m . 




m (g) 

Fig. 1. — The mass dependence, $(m), of the adopted model for the collisional lifetime of 
particles. The plot shows $(m) for a particle on the circular orbit at the reference heliocentric 
distance R 0 = 1 AU. Two parameters of the collisional model, Si and S 2 , were used to test 
the sensitivity of our results to modifications of $(m). 
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Fig. 2. — Effects of P-R drag on the highly-eccentric orbits of OCC particles, (a) Evolution 
tracks of particles evolving by P-R drag in semimajor axis and perihelion distance. Particles 
with e ~ 1 evolve from right to left along the lines of constant perihelion distance. The 
shaded area is inaccessible to orbits, (b) Evolution timescale for particles with D = 100 /mi 
and p = 2 g cm -3 . Contours show the time of fall, Tf a n, from the initial orbit with a and q 
to the Sun. According to Eq. (11), Tf a n scales linearly with D (and p), so that, for example, 
the 1 My contour for D = 100 pm is the 10 My contour for D = 1 mm. 
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Fig. 3. — The orbit history of a particle with D = 100 fin i, initial a = 10 3 AU and q = 
1.53 AU. After 1.7 Myr, the orbit decouples from Jupiter and moves to a ~ 1 AU. The 
particle ends up having a < 0.1 AU, and sublimates upon reaching R < 0.05 AU. 
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Fig. 4. — The radial distribution of OCC particles that we obtained in our model. The 
number density was normalized to the input flux of one particle released on bound orbit 
per year. We used D = 100 /mi here. The distribution with D = 300 /irn looks similar 
for R > 5 AU, but is depleted below 5 AU, relative to the one shown here, because fewer 
particles with D = 300 /irn are able to decouple from Jupiter. Particles were released on 
orbits with a ~ 10 3 AU and uniformly random sin i. The initial distribution dN(q) was set 
as described in Section 2.1 with 7 = 0.5 (the results are not sensitive to 7 ). 
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Fig. 5. — The steady state distribution of orbital elements of particles with D = 100 jin i: 
(a) d N(a), (b) diV(e), and (c) d N(i). The distributions shown here ignore particles with 
R > 5 AU. They therefore represent the steady state distribution of orbits in the inner solar 
system, which is relevant for observations of the inner zodiacal clond and sporadic meteors. 
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Longitude (deg) 


Fig. 6. — The radiants of OCC meteoroids determined in our model with I* = 0.01: (a) 
D = 100 jum, and (b) D = 300 /mi. The north and south apex sources are clearly visible. In 
addition, the radiant distribution in panel (b) shows a ring structure centered at (/, b) = (0, 0) 
and extending to ~60° in l and b. The high-frequency fluctuation of radiant density between 
neighbor bins is due to insufficient statistics and should be disregarded. The units of the 
side bar are arbitrary. 
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Fig. 7. — The same as Fig. 6 but for D = 100-/mr particles released from lP/Halley. The 
results for 55P/Tempel- Tuttle are similar. 
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Fig. 8. — The orbital element distribution of model apex meteoroids: (a) v g , (b) a, (c) i, and 
(d) e. Here we used OCC particles with D = 100 /mi and I* = 0.01. The apex meteoroids 
were selected by using the following radiant cutoffs: —40° < l < 40° and —40° < b < 40°. 
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Fig. 10. — The expected relative strength of the apex and helion/antihelion sources as a 
function of particle size. The triangles (diamonds) show the result for OCC (lP/Halley) 
particles. We calculated the relative strength by dividing Pi(D) listed in Table 1 for OCC 
(lP/Halley) particles (assumed here to represent the apex source) by Pi(D) of 2P/Encke 
particles (assumed to represent the helion/antihelion source; W09). The Pi(D ) ratio is a 
proxy for the relative strength of meteor sources if dNo(D) of different initial populations 
had roughly the same shape. We arbitrarily normalized the ratio to 1 for particles with 
D = 100 /jrn. As more sensitive radars detect smaller particles, they are expected see more 
apex meteors, because Pi(D ) of OCC (lP/Halley) particles increases while that of JFC 
particles drops. 




